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^ . Abstract 

Quantum field theories in front-form dynamics are not manifestly rotationally in- 
, variant. We study a model bound-state equation in 3+1 dimensional front-form 

dynamics, which was shown earlier to reproduce the Bohr and hyperfine structure 
of positronium. We test this model with regard to its rotational symmetry and find 
^ . that rotational invariance is preserved to a high degree. Also, we find and quantify 

the expected dependence on the cut-off. 

o 
in 
o 



^ ■ The framework of Discretized Light Cone Quantization (DLCQ) has been 

applied successfully to many and diverse physical systems, cf. Ref. [1] and 
references therein. The method is especially efficient when applied in lower 
Q^l dimensions, and seems ideally suited for strictly two-dimensional systems [2,3]. 

The extension of the DLCQ program into the physical four dimensions is both 
necessary and exciting, even when some of its striking advantages are lost on 
the way. Several attempts to do so were hampered by a prominent draw-back 
rS I of the front form, particularly its lack of manifest rotational symmetry; but 

I recently some progress has been made [4]. 

The method of DLCQ was applied to 3+1 dimensions first by Tang et al. [5], 
to check the method at the example of positronium. This approach, like others 
[6] trying to solve the matrix eigenvalue problem by using a light-cone adapted 
Tamm-Dancoff procedure, suffered severe convergence problems. Wilson and 
collaborators [7] considered the problem in more abstract terms, emphasiz- 
ing the role of renormalization problems. But one may state in all fairness 
that concrete and practical prescriptions have not emanated thus far. Their 
positronium spectrum [8] was obtained later. Within the model of Krautgart- 
ner, Pauli and Wolz [9], however, the Bohr and the hyperfine structure of 
positronium on the light-front was resolved, but only for the z-component of 
total angular momentum = 0. The derivation of the effective interaction 
was not overwhelmingly convincing, and the problems with a cutoff-dependent 
eigenvalue equation have not been settled thus far, see however also Ref. [10]. 
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But even with the apparent shortcomings of this model, one can examine 
its rotational symmetry by asking whether the corresponding members of a 
rotational multiplet are degenerate or not, and if the multiplets contain the 
correct number of degenerate states. The agenda is then quite obvious: 

• Generalize the method to arbitrary J^; 

• Increase the numerical accuracy needed for that. 

If we can calculate the eigenvalue spectrum separately for each J^, which is 
a kinematic operator both in the front and in the instant form [11], some 
of the eigenvalues of distinct must be degenerate and form a multiplet, 
also in the front form. We can thus investigate quantitatively to which extent 
rotational symmetry is violated within a specific model. Any violation of the 
degeneracy of multiplets is a strong indication that the original covariance of 
the Lagrangian was lost by nature of the (model dependent) approximations. 



1 Model positronium 

The model of positronium considered here has been introduced by Krautgart- 
ner et al. [9] and we refer to it for all unquoted details. In light-cone quantiza- 
tion, the contraction of the momentum operators P'^, is called the light-cone 
Hamiltonian, Hi^c — P^Pn- Solving the eigenvalue problem 



yields the mass (squared) eigenvalue spectrum of a physical system. This full 
problem is very difficult to solve, particularly in gauge theory, see [12]. It is 
easier to solve the problem in a reduced space with an effective Hamiltonian 



with some free kinetic part T and an effective interaction U . The reduced 
space is here the Fock space of the single electron and a single positron, i.e. 
1'^) = l^'ee)- In what follows, we denote the mass, the longitudinal momen- 
tum fraction, the transverse momentum, and the helicity of the electron with 
m, X, A?L, and Ag, respectively, and those of the positron with m, 1 — x, —k±, 
and Ag. Physically, the effective interaction scatters an electron-positron pair 
from a state with four-momenta {ke,ke) into a state with {k'^,kg), which in 
general has a different free invariant mass than the entrance channel. This 
requires a certain amount of symmetrization which was made in all previous 
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work [14,15], and thus also in [9]. A more thorough discussion can be found 
Ref. [17]. This in mind, Eq. (2) becomes an integral equation 



M^{x,kr,XeAe\'ipn) 




{x, k±; Ae, 



(3) 



x{l — x) 



27r2 



a 




dx'd?k'^ 



{x, kr, Ae, X,\S\x', fci; a;, X',){x\ k[; A^, A;|^„). 



Symmetrization is reflected in the denominator 



Q"^ {x , k±; x , kj_) — — {ke k^) + {ke k^) , 



(4) 



which is the mean Feynman four-momentum transfer. The Lorentz-contracted 
Dirac spinors are collected in 



(Ae, Ae-|5| A'e, A^-) = [u(k,, K)l'u{k',, A'e)] [v{k',, X'^^Ak^, A.)] ■ (5) 



In hehcity space it is a 4 x 4 matrix which is tabulated exphcitly in the 
Compendium [16] or in Ref. [18, App. B]. The domain of integration D is 
specified by the 'sharp' cut-off A 



i.e. by Lepage-Brodsky regularization [15]. The experimental literature differ- 
entiates between proper positronium (//e) and true positronium (ee) (Telegdi). 
The model in Eq.(3) is neither of the two: It is 'model positronium'. The an- 
nihilation channel is however included in Refs. [17-19]. 

In the meantime, the inherent difficulties with the light-cone adapted Tamm- 
Dancoff approach have faded away, since better formalisms yield the very same 
integral equation (3), see f.e. [10]. Note that the above equation is actually a 
set of four coupled integral equations in the three variables x, k^ — \k± \ simp 
and ky = \k± \ cos ip. Its numerical solution is highly non-trivial. 

2 The construction of the multiplets and the numerical results 

From usual (instant form) quantization one is used to the fact that the Hamil- 
tonian commutes with the total angular momentum operator with eigenval- 
ues J(J + 1). As a consequence, the Hamiltonian eigenvalues are (2J + l)-fold 



+ k[^ < x'{l - x')(A2 + 4m2), 
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Fig. 1. The spectrum of the positronium model for | J^l < 3 at A = m, Ni=N2=21. 

The eigenvalues M are given in units of the electron mass m. - The multiplct struc- 
ture of the eigenvalues is emphasized by the phenomenological notation '^^^^Lj". 



degenerate. The individual members of a J^-multiplet are labeled by J^. In the 
front form things are different: Jz is the only kinematic of the three rotation 
operators. The key observation here is that front-form eigenstates can be clas- 
sified by Jz but not by a quantum number with respect to J^. We can, however, 
diagonalize the Hamiltonian subsequently in different sectors characterized by 
Jz, i.e. as a function of Jz = 0, ±1, ±2, . . .. Some of the eigenvalues will be 
degenerate in Jz, within a certain numerical accuracy of course, and form the 
usual (instant form) multiplets. The largest value of in a given multiplet 
determines the multiplicity, 2{Jz)mayi + 1; and thus indirectly the quantum 
number J. If we find these degenerate multiplets, we have reproduced the 
multiplet structure of the instant form. The question is then to which extent 
multiplets are degenerate, i.e. on which level of accuracy rotational invariance 
is violated. 
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Our prejudices arc illustrated in Fig. 1. Here is our most important result, and 
we shall explain below how to get it. The figure has the anticipated proper- 
ties: the eigenvalues are arranged in degenerate multiplets, and each multiplet 
has an odd number of members. The eigenvalues can be arranged in clusters 
which are characterized by the Bohr quantum number n. For a pure Coul- 
omb spectrum there would be (4n — 2) degenerate states with spin projection 
Jz = 0, and we found that the multiplets in our calculations pass this check: 
they have precisely (4n — 2) members. We have observed this property up to 
n—5, beyond which we have seen no reason to pursue. 

Another important result is that only those combinations of the quantum 
numbers ttc and vr^^ appear which are expected from a non-relativistic analysis. 
Usually one classifies the states with the spectroscopic notation ^'^"'"^Lj. In 
front-form dynamics, however, neither the total angular momentum J, nor 
the orbital angular momentum L, nor the total spin S are good (kinematic) 
quantum numbers. Rather, the front-form Hamiltonian is symmetric under 
charge conjugation C and a combination of time reversal T and parity V, 
sometimes called handedness 7i = VT . In the non-relativistic case, one can 
relate the (instant form) quantum numbers (J, L, S) uniquely to the (front 
form) quantum numbers (ttctt-^), if one chooses a certain convention for the 
time reversal operation [20]. The spectroscopic notation is used only to label 
the states conveniently [17-19]. 

The multiplets of mass (squared) eigenvalues are seemingly degenerate on the 
scale of the figure (Fig. 1). To be more quantitative we collect in Table 1 the 
first 18 eigenvalues in the form of binding coefficient defined by 

B^{J,)^^ {2m -Mr,), (7) 

where the trivial dependence on mass and coupling constant is removed. The 
numerical errors (in parenthesis) are estimated from the difference between 
the values for maximal and next to maximal number of integration points. 
Except for a few states to be discussed below, the relative discrepancy of 
corresponding eigenvalues is typically a few parts in 10^, i.e. the degeneracy 
is on the level of numerical errors of the diagonalization routine. 

The positronium spectrum has been calculated perturbatively, long ago by 
[21] and [22], and we have to keep in mind that these calculations were done 
under the proviso of the small a ~ 1/137. Despite that, our eigenvalues agree 
quite well, i.e. on the percent level, with those results, even for the highly 
excited states, see Table 1. 

The singlets tend to be too weakly bound, especially for the low Bohr 
quantum numbers n — 1 and n — 2. This effect is reversed for the triplets ^Si. 
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Table 1 

The binding coefficients Bn for a = 0.3, A = m, A^i = N2 = 21. — Column 1 
gives the spectroscopic notation n ^^~^^Lj and column 2 the parities C and H. The 
numerical errors of 5„ (in parenthesis) are estimated, sec text. The last column lists 
the discrepancy to perturbation theory up to order 0{a'^) [22]. 
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0.25526(18) 
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0.11344(36) 
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0.11298(16) 


0.11298(16) 


-0.43 
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-1.03 



Also, the ordering of the multiplets seems to have minor errors. For instance, 
the 2 ^5*0 state and the 2 ^Pq state are permuted: the -S-state should be the 
lowest according to perturbation theory. We convinced ourselves that this is 
a finite cut-off effect. 



The numerical approach has two formal parameters: the number of integration 
points Ni — N2 = N and A. The dependence on the number of integration 
points N was found to decrease exponentially fast, even for the splitting be- 
tween triplet states. The asymptotic value AM^(A^^oo) = a amounts to 
only 0.5% of the relevant scale, namely the singlet-triplet splitting, of roughly 
0.0102m^. With other words, the states are highly degenerate. 
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Table 2 



The binding coefficients of the singlet (Bg) and the triplet states (Bt) for a = 0.3, 
A'^i = 25, N2 = 21 are given as function of the cut-off A. They are compared with 

the r(^sulls of p(MMur})ali()ii lli(M)ry [22] uj) lo order 0{c\ ^) ai:d ©(n'Mrin). 



-J^ 

cut-off — 
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Bt 




1.0 
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10.8 


1.37112216 
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1.38744792 
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16.2 
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18.0 


1.42740143 


1.13641774 


0.80828803 


0{a^ In a) 


1.11812500 


0.99812500 


0.33333333 
0.23792985 



The dependence on the cut-off A. The weak point of the present work 
is the expected dependence on A. Its occurrence is not a particular feature 
of the light cone approach, but appears due to the Dirac interaction in any 
(gauge) field theory: For very large (transversal) momenta, i.e. for k'^ ^ k^, 
holds (ti \S\ U)/Q^ — ' 2 and (it \S\ iD/Q^ — > 2: The kernel does not 
decay sufficiently fast, see also Ref. [10]. 

For discussing the dependence on A to some detail wc present in Table 2 the 
binding coefficient for the singlet and the triplet together with the hyperfine 
coefficient 

^hf = J {Mtriplet — Msinglet) ■ (8) 



Its perturbative value [21,22] is known: 



^('"2 + y)+o(.M„„). 



(9) 



The term in square brackets is the contribution from the one-photon anni- 
hilation. When comparing to previous results [9, Table V], we find that the 
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Fig. 2. The discrepancy function ^Mf 
(= M2(1) - M2(0)), is plotted versus the 
cut-off A for the V^Si (A) [upper plot], 
the 2^Pi, and the 2^Pi (•) [lower plot]. 
Parameters are A = m, A/^i = 21, and 
N2 = 21. 



Fig. 3. The decrease of the = 1 triplet 
ground state wave function with paral- 
lel helicities function of the off-shell 
mass //. The parameters are A = 20.0 m, 
A^i=41, A^2=ll- The six different curves 
correspond to six values of the angle 6. 
Rotational symmetry is only for small /i. 



singlet falls off faster with A, yielding a smaller value for the hyperfine coeffi- 
cient in our calculations. We note that the predictions of perturbation theory 
for different orders differ as much as 40% for the value of a = 0.3 used in the 
present work. 

The degeneracy of eigenvalues as a function of the cut-off A is displayed in 
Fig. 2. One notices that the discrepancy between the triplet ground states 
rises from close to zero at A = m to roughly 50% of the hyperfine splitting for 
the large value of A = 50m. For the excited states the difference stays always 
below 12% on the relevant scale. Since we expect a logarithmic divergence 
we fit the curves M^(A) to a polynomial in In A. If one omits the points for 
A > 20 m because the integrations of the numerical counter term can become 
problematic in this region, the fit yields 



Mlnaieti^) = i 3-90545 - 0.03510 In — + 0.00746 In^ —\ m^, 
" \ m mJ 

Mtripieti^) = (^3.90976 - 0.01858 In — + 0.00789 In^ — ) 

The small coefficient of the In^A/m-term hints thus indeed at a logarithmic 
cut-off dependence of the eigenvalues. The dependence weakens if one includes 
the annihilation channel [19]. 



(10) 



3 More details for general Jz and the wavefunction 

We want to calculate the Hamiltonian spectrum in all sectors of J^. It is ad- 
vantageous to reformulate the eigenvalue problem in such a way that allows 
for a separation of the sectors, i.e. a block diagonalization of the Hamilto- 
nian with respect to the conserved quantum number Jz. By exploiting the 
symmetry of the Lagrangian under rotations in the (x, y) plane, we can sub- 
stitute the angular variable </? by the quantum number Jz applying a integral 
transformation on Eq.(3), and in particular on the effective matrix elements 

(X, k^, Jz-, Xe, XelU'^'lx', fcl, Jz-, K, K) (11) 
27r 2n 

= —J jd^d^' e-^(^^^-^i'^')(x, k^, ^- Ae, Ae-|C/^^|x', K, X',). 



It is easy to convince one-self that angular dependence can enter only in the 
form of the difference cp — ip', and that the matrix elements are well-behaved 
for an arbitrary Lz—Jz—Sz- The eigenvalue problem, Eq. (3), looks now like 

0=(^M'n- Ae, Ae-; Jz^ (12) 

+ dx'dk'^{x, k^- Ae, X,\U^^\x\ fcj; X'^, X',) {x', fcj; A^, A'^; J,|^„). 

The angle averaged 'general helicity table' (Ag, Xe\U'^^\X'^, Ag) is again a 4 x 4- 
matrix in helicity space and tabulated in Refs. [17-19]. 

In the pioneering work of Kraut gartner et al. [9] Eq. (12) was solved numer- 
ically for the special case = 0. By convenience they use momentum coor- 
dinates {/J,, 9) instead of {x,k_i), corresponding to the Sawicki transformation 
in the Compendium [16]. The integral equation (3) is converted into a matrix 
equation by Gauss-Legendre quadratures in the off-shell mass n — 2\k\ and 
the polar angle 9, with A^i and A^2 integration points, respectively. The Hamil- 
tonian matrix is subsequently diagonalizcd by numerical methods. In Rcf. [9] 
the fine structure constant was set to the very large value a = 0.3 in order to 
resolve the accumulation of the eigenvalues around ~ (2m) ^. We use the 
same unphysically large coupling to facilitate the search for possible violations 
of rotational symmetry in the spectrum. Because of the inherent quadratically 
intcgrable Coulomb singularity (Q"^) convergence of the eigenvalues with the 
number of the (Gaussian) integration points is extremely slow. Crucial im- 
provement is achieved by using the Nystr0m method [23]. Its essence is to add 
a diagonal numerical (Coulomb) counter term and to subtract its discretized 
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Fig. 4. Triplet wavefunction ^TT ^ -^i§- 5- Triplet wavefunction ■i/'ij is for 
function of x and |A-_l|, at = 0. = 1, A = m, iVi = 41, and iV2 = 11. 

version. In contrast to the pure Coulomb case, where the only counter term can 
be calculated analytically [24], the fine and hyperfine interactions in Eq.(3) 
require several counter terms. For J^=0 different diagonal matrix elements 
occur for parallel and anti-parallel helicities; they require two counter terms. 
Since both amplitudes have a similar singularity structure and comparable 
values, Krautgartner et al. have used the same semi-analytical counter term 
in both diagonal matrix elements. Despite of this unnecessary simplification 
the convergence of the eigenvalues was satisfactory. In order to analyze the 
multiplet structure of the spectrum in all sectors of J^, we have to improve 
the accuracy of this numerical algorithm. When J 7^ 0, one has four different 
diagonal elements in the helicity matrix, one of which is much smaller than the 
others. To meet these numerical requirements, we avoid analytical treatment 
of the counter terms altogether, and integrate the counter terms numerically 
with a sufficiently high precision. This requires a somewhat larger computing 
time, but one is rewarded by an improvement of a factor two in convergence. 

The wavefunctions are generated simultaneously with the spectrum. The two 
components of the lowest state [J^ = 0) with anti-parallel and parallel he- 
licities, look, up to a scaling factor, like the wavefunctions with parallel and 
anti-parallel helicities of Figs. 4-7, respectively. The shapes and the peak values 
are the same as in Ref. [9]. The plots in Ref. [9] seemed to indicate numerical 
problems because they showed internal structure. We found that this is due 
to numerical mistakes in the graphing package. 

Let us discuss briefly the properties of the Jz=^ wavefunctions, as displayed 
in Figs. 4-7. Due to the lower symmetry, the wave functions for 7^ show 
more structure than those for = 0. The wave functions with 7^ have 
four components corresponding to the four different helicity combinations. We 
can see immediately from Figs. 5 and 6 that the components for anti-parallel 
helicities are identical. The components for parallel helicities in Figs. 4 and 
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Fig. 6. Triplet wavefunction ip^p 



Fig. 7. Triplet wavefunction -0^. 



7 have rather disjunct properties: the (tt)-component peaks at a; = 0.5 and 
k± = and is rotationally invariant. The (||)-component vanishes at = 
and is shaped more like the components with anti-parallel helicities. Note 
the very different peak values: the anti-parallel components are suppressed 
by a factor of 40 as compared to the (tt)-component, the (||)-component is 
suppressed by a factor of 1400! We emphasize that the wave functions are 
not rotationally invariant, despite the fact that their non-relativistic analogue 
is an s-wave. One reason is the fact that rotations around the transverse 
axes are dynamical. Another way of seeing this is to transform the light-cone 
variables to quasi-equal-time coordinates. The Jacobian of the transformation 
breaks rotational invariance due to the term kl / {rn^ + k1 + k'^) . The breaking of 
rotational invariance of the wave function is noticeable only for large momenta 
^ m, or, correspondingly, large cut-offs A. As an example, the triplet wave 
function l'^S'i(tt) is plotted in Fig. 3 as a function of the off-shell mass /x. 
Obviously, the wave function is not isotropic but depends on the angle 6 for 
sufficiently large values of the off-shell-mass. We note that with a cut-off of 
A = m in place, the discrepancy of curves with different cos^ would not be 
resolvable. 



Conclusions: 

(1) The counter term technology of Krautgartner el al. [9] can be implemented 
without unnecessary simplifications. 

(2) The generalization to all values of the angular momentum projection is 
possible and novel. 

(3) The dependence on the cut-off is logarithmic, and thus renormalizable. 

(4) Possible problems with rotational symmetry can be dealt with easier in 
practice than anticipated by more formal investigations in the literature. 

(5) The eigenvalues for different Jz arrange themselves in highly degener- 
ate multiplets. This in turn lets us conclude that the effective interaction of 
Krautgartner et al. [9] is consistent with the rotational thus Lorentz-invariance 
of the Lagrangian, and this holds probably for all light-cone based approaches. 
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